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We study the heat kernel for a Laplace type partial differential operator 
acting on smooth sections of a complex vector bundle with the structure 
group GxU(\) over a Riemannian manifold M without boundary. The total 
connection on the vector bundle naturally splits into a G-connection and a 
i7(l)-connection, which is assumed to have a parallel curvature F. We find 
a new local short time asymptotic expansion of the off-diagonal heat kernel 
U(t\x, x') close to the diagonal of M x M assuming the curvature F to be 
of order t~^. The coefficients of this expansion are polynomial functions in 
the Riemann curvature tensor (and the curvature of the G-connection) and 
its derivatives with universal coefficients depending in a non-polynomial 
but analytic way on the curvature F, more precisely, on tF. These func- 
tions generate all terms quadratic and linear in the Riemann curvature and 
of arbitrary order in F in the usual heat kernel coefficients. In that sense, 
we effectively sum up the usual short time heat kernel asymptotic expan- 
sion to all orders of the curvature F. We compute the first three coefficients 
(both diagonal and off-diagonal) of this new asymptotic expansion. 
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1 Introduction 

The heat kernel is one of the most powerful tools in quantum field theory and 
quantum gravity as well as mathematical physics and differential geometry (see 
for example [[lilllSlElinillSllISlEZllriland further references therein). It is 
of particular importance because the heat kernel methods give a framework for 
manifestly covariant calculation of a wide range of relevant quantities in quantum 
field theory like one-loop effective action, Green's functions, effective potential 
etc. 

Unfortunately the exact computation of the heat kernel can be carried out 
only for exceptional highly symmetric cases when the spectrum of the operator 
is known exactly, (see [[I71[l9l|20] and the references in [?l[l5l[l4l[l3l). Although 
these special cases are very important, in quantum field theory we need the ef- 
fective action, and, therefore, the heat kernel for general background fields. For 
this reason various approximation schemes have been developed. One of the old- 
est methods is the Minackshisundaram-Pleijel short-time asymptotic expansion of 
the heat kernel as ? ^ (see the references in lfT8l l2ll23ll'). 

Despite its enormous importance, this method is essentially perturbative. It 
is an expansion in powers of the curvatures R and their derivatives and, hence, is 
inadequate for large curvatures when tR ~ 1. To be able to describe the situa- 
tion when at least some of the curvatures are large one needs an essentially non- 
perturbative approach, which effectively sums up in the short time asymptotic 
expansion of the heat kernel an infinite series of terms of certain structure that 
contain large curvatures (for a detailed analysis see [4, 9| and reviews [lOlfTllI '). 
For example, the partial summation of higher derivatives enables one to obtain 
a non-local expansion of the heat kernel in powers of curvatures (high-energy 
approximation in physical terminology). This is still an essentially perturbative 
approach since the curvatures (but not their derivatives) are assumed to be small 
and one expands in powers of curvatures. 

On another hand to study the situation when curvatures (but not their deriva- 
tives) are large (low energy approximation) one needs an essentially non-pertur- 
bative approach. A promising approach to the calculation of the low-energy heat 
kernel expansion was developed in non-Abelian gauge theories and quantum grav- 
ity in [i3i|4i [51 [6l El [Hill [Hill. While the papers USUI [6113 dealt with the paral- 
lel f/(l)-curvature (that is, constant electromagnetic field) in flat space, the papers 
im [8l [131 dealt with symmetric spaces (pure gravitational field in absence of an 
electromagnetic field). The difficulty of combining the gauge fields and gravity 
was finally overcome in the papers [[T4l [TSl . where homogeneous bundles with 
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parallel curvature on symmetric spaces was studied. 

In this paper we compute the heat kernel for the covariant Laplacian with a 
large parallel U{1) curvature F in a Riemannian manifold (that is, strong covari- 
antly constant electromagnetic field in an arbitrary gravitational field). Our aim is 
to evaluate the first three coefficients of the heat kernel asymptotic expansion in 
powers of Riemann curvature R but in all orders of the C/(l) curvature F. This is 
equivalent to a partial summation in the heat kernel asymptotic expansion as f — > 
of all powers of F in terms which are linear and quadratic in Riemann curvature 
R. 

2 Setup of the Problem 

Let M be a n-dimensional compact Riemannian manifold without boundary and 

be a complex vector bundle over M realizing a representation of the group 
G®U{1). Let ^ be a section of the bundle S and V be the total connection on the 
bundle S (including the G-connection as well as the U(l)-connection). Then the 
commutator of covariant derivatives defines the curvatures 

[V^,V,]^ = (!??^, + /F^>, (2.1) 

where !R^y is the curvature of the G-connection and F^y is the curvature of the 
?7(l)-connection (which will be also called the electromagnetic field). 

In the present paper we consider a second-order Laplace type partial differen- 
tial operator, 

^ = -A, A = /^V,Vy . (2.2) 
The heat kernel for the operator ^ is defined as the solution of the heat equation 

(d, + ^)U(t\x,x') = , (2.3) 

with the initial condition 

U(0\x, x') = P(x, x')6(x, x') . (2.4) 

where S(x, x') is the covariant scalar delta function and P(x, x') is the operator of 
parallel transport of the sections of the bundle S along the geodesic from the point 
x' to the point x. 

The spectral properties of the operator =Sf are described in terms of the spectral 
functions, defined in terms of the traces of some functions of the operator =Sf , 
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such as the zeta-f unction ^{s) = Tr ^ \ and the heat trace 
Tr expC-?^) = ^ Jvol tr W^'^^it) , 



(2.5) 



M 



where Jvol = g^^^dx is the Riemannian volume element with g = det^^y and tr 
denotes the fiber trace. Here and everywhere below the diagonal value of any two 
point quantity /(jc, x') denotes the coincidence limit as x ^ x' , that is, 

/-g = /(^,x). (2.6) 

It is well known IfTSl that the heat kernel has the asymptotic expansion as 
t (see also |I21[IUI[II1|311) 

U{t\x, x) ~ {AnlV'^nx, x')h}'^{x, x) exp 



(t{x, x') 

2r~ 



J^t'atix,x') , (2.7) 



k=0 



where cr(x, x') is the geodesic interval (or the world function) defined as one half 
the square of the geodesic distance between the points x and x' and A{x, x') is 
the Van Vleck-Morette determinant. The coefficients a^ix, x') are called the ofF- 
diagonal heat kernel coefficients. 

The heat kernel diagonal and the heat trace have the asymptotic expansion as 



k=0 



where 



and 



k=0 

00 

Tr exp(-?^) ~ 

diag f N 

<^k ^k{X, X) 

Ak = d\ol tr af 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



M 



The coefficients A^ are called the global heat kernel coefficients; they are spectral 
invariants of the operator =Sf . 
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The diagonal heat kernel coefficients a^'^^ are polynomials in the jets of the 
metric, the G- connection and the i7(l) -connection; in other words, in the curva- 
ture tensors and their derivatives. Let us symbolically denote the jets of the metric 
and the G-connection by 

Rin, = |V,^. ■ • • V,„7?V.> v..) . V(,. • ■ ■ V,X,.,.o] . (2.12) 
and the jets of the i7(l) connection by 

i^(n) = V(„ ■ ■ ■ V, FVo ■ (2.13) 

Here and everywhere below the parenthesis indicate complete symmetrization 
over all indices included. 

By counting the dimension it is easy to describe the general structure of the 
coefficients af'^^. Let us introduce the multi-indices of nonnegative integers 

i = (zi,...,zj, j = (2.14) 

Let us also denote 

|i| = + . . . + = + . . . + . (2.15) 

Then symbolically 



a, 



_k N N-l 

diag 



^ ^ ^ ^ C(k,l,m)XiP'(ji) ■ • ■ -^"(7;) -^('i) ■ ■ ■ R(i„,) ' (2-16) 



^=1 /=0 ;;i=0 iJ^O 

|l|+lil+2«=2*' 



where C(k,i,m)X} some universal constants. 

The lower order diagonal heat kernel coefficients are well known [|T8l l2l[Tn 

= 1, (2.17) 



af'' = ^R, (2.18) 



2 30 72 180 180 

+^'R,yW' + I'R.JF^'' - ^F^^F^' . (2.19) 



To avoid confusion we should stress that the normalization of the coefficients aj, 
differs from the papers El fTOlfTTll . 
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In the present paper we study the case of a parallel U{\) curvature (covariantly 
constant electromagnetic field), i.e. 

V^F,/^ = . (2.20) 

That is, all jets F(„) are set to zero except the one of order zero, which is F itself. 
In this case eq. (12.161) takes the form 

*: N N-l 

^ ^(*^''.'").''^' -^('l) ■ "^iim) ' (2.21) 

N=l 1=0 m=0 iao 

li|+2A'=2* 

where C^k,i,m),\ are now some (other) numerical coefficients. 

Thus, by summing up all powers of F in the asymptotic expansion of the heat 
kernel diagonal we obtain a new (non-perturbative) asymptotic expansion 

where the coefficients a^'^^(0 are polynomials in the jets i?(„) 

<"« = ZZ Z •••%.)' (2-23) 



^=1 m=0 i>0 

|i|+2«=2/S 



and fl^i-ft) are some universal dimensionless tensor-valued analytic functions that 
depend on F only in the dimensionless combination tF. 

For the heat trace we obtain then a new asymptotic expansion of the form 

DO 

Tr exp(-?^) ~ {AntV'^ Yj ^^k{t) , (2.24) 

k=Qi 

where 



^ dvo\ tra^'^^ 

M 



Au{t) = dvol tra,\t) . (2.25) 



This expansion can be described more rigorously as follows. We rescale the 
i7(l) -curvature F by 

F ^ F{t) = t-^F , {126) 
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so that tF{t) = F is independent of t. Then the operator J^(t) becomes dependent 
on t (in a singular way!). However, the heat trace still has a nice asymptotic 
expansion as t ^ 

CO 

Tr exp[-t^(t)] ~ {4ntr"'^ , (2.27) 

where the coefficients Aj, are expressed in terms of F = tF{t), and, therefore, are 
independent of t. Thus, what we are doing is the asymptotic expansion of the heat 
trace for a particular case of a singular (as t ^ 0) time-dependent operator =Sf (f). 

Let us stress once again that the eq. (12.231) should not be taken literally; it only 
represents the general structure of the coefficients af'^^it). To avoid confusion we 
list below the general structure of the low-order coefficients in more detail 

at\t) = fyt), (2.28) 
af^'it) = /(Vj;^'^^(0^./?,v + /S)''''(0^pv, (2.29) 

, r(2) a/iySnva-p f .\n p . f(2) afitiV(.\fn (n 

■'(2,1) \'' J^afSyS^fivcrp J (2,2) KUJ^afS'^iiv 

+/g;^^-''(OV.v.p (2.30) 

with obvious enumeration of the functions. It is the universal tensor functions 
fg]jj^(t) that are of prime interest in this paper. Our main goal is to compute the 
functions fgl,)(t) for the coefficients afj' ^^(0, af'^^{t) and df'^^it). 

Of course, for ? = (or F = 0) the coefficients ak(t) are equal to the usual 
diagonal heat kernel coefficients 

af'\0) = af'''. (2.31) 

Therefore, by using the explicit form of the coefficients af'^^ given by (12.191 ) we 
obtain the initial values for the functions fj^j^j^y Moreover, by analyzing the corre- 
sponding terms in the coefficients af'^^ and a'^'^^ (which are known, [|T8l |2l l22l). 
one can obtain partial information about some lower order Taylor coefficients of 
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the functions f^jk)(f)'- 



fyt) = l-lt%yF^' + 0(t'), (2.32) 
/(iW(0 = ^5'^<] + O(0, (2.33) 
/(lirW = ItiF^' + Oif), (2.34) 



fSi"'''^<rp(t) = lg«^5^5;j + O(0, (2.35) 



/(?2;V(o = -^tiF^^^ + oit"), 



(2.36) 



J{2,\)«P l^y j^gQ<5^[a<5/3]vS S j^gQ^ [0-5/3] [v6 

(2.37) 

/(S;V(0 = j^^<j + O(0, (2.38) 

(2.39) 

This information can be used to check our final results. 

Notice that the global coefficients A^it) have exactly the same form as the local 
ones; the only difference is that the terms with the derivatives of the Riemann cur- 
vature do not contribute to the integrated coefficients since they can be eliminated 
by integrating by parts and taking into account that F is covariantly constant. 

Moreover, we study even more general non-perturbative asymptotic expansion 
for the off-diagonal heat kernel and compute the coefficients of zero, first and 
second order in the Riemann curvature. We will show that there is a new non- 
perturbative asymptotic expansion of the off'-diagonal heat kernel as ? ^ (and 
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F = t ^F, so that tF is fixed) of the form 



U(t\x, x) ~ r{x, x)^^'^{x, x)UQ{t\x, x) ^ t^'^bkiAx, x) (2.40) 



k=0 



where f/o is an analytic function of F such that for F = 

o-{x, x') 



UQ{t\x,x') 



= {Antyl^ exp 



F=0 



2t 



(2.41) 



Here bk{t\x, x') are analytic functions of t that depend on F only in the dimension- 
less combination tF. Of course, for ? = they are equal to the usual heat kernel 
coefficients, that is. 



b2k{0\x, x) = akix, x) , b2k+\{0\x, x) = 0. 



(2.42) 



Moreover, we will show below that the odd-order coefficients vanish not only for 
t = and any x + x' but also for any t and x = x\ that is, on the diagonal, 

4!fj(0 = 0. (2.43) 

Thus, the heat kernel diagonal has the asymptotic expansion (|2.22l) as ? — > with 

~af\t) = {Antr''ui''''{t)bf^\t) . (2.44) 



3 Geometric Framework 

Our goal is to study the heat kernel U{t\x, x') in the neighborhood of the diagonal 
as X ^ x'. Therefore, we will expand all relevant quantities in covariant Taylor 
series near the diagonal following the methods developed in |12l[TTl[T0l[l2l|. We fix 
a point, say x', on the manifold M and consider a sufficiently small neighborhood 
of x', say a geodesic ball with a radius smaller than the injectivity radius of the 
manifold. Then, there exists a unique geodesic that connects every point x to 
the point x'. In order to avoid a cumbersome notation, we will denote by Latin 
letters tensor indices associated to the point x and by Greek letters tensor indices 
associated to the point x'. Of course, the indices associated with the point x (resp. 
x') are raised and lowered with the metric at x (resp. x'). Also, we will denote by 
Vq (resp. Vp covariant derivative with respect to x (resp. x'). 
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We remind below the definition of some of the two-point functions that we 
will need in our analysis. First of all, the world function cr(x, x') is defined as one 
half of the square of the length of the geodesic between the points x and x'. It 
satisfies the equation 

cr = -u^Ua = u^u^ , (3.1) 

where 

Ua = ^aO-, Uf, = V^Cr . (3.2) 

The variables are nothing but the normal coordinates at the point x'. 
The Van Vleck-Morette determinant is defined by 

A(x,x') = g-Hx)dct[-VaY,o-(x,x')]g-kx') . (3.3) 

This quantity should not be confused with the Laplacian A = g'^^V^Vy. Usually, 
the meaning of A will be clear from the context. We find it convenient to parame- 
terize it by 

A(x, x') = exp[2^(x, x')] . (3.4) 

Next, we define the tensor 

rfh = V.Y'^o- , (3.5) 
and the tensor 7"^ inverse to rfa by 

y%rfb = 6l, rf,y\ = ^,. (3.6) 

This enables us to define new derivative operators by 

V;. = y%^a . (3.7) 

These operators commute when acting on objects that have been parallel trans- 
ported to the point x' (in other words the objects that do not have Latin indices). 
In fact, when acting on such objects these operators are just partial derivatives 
with respect to normal coordinate u 



We also define the operators 



1)/. = - :^iF,au" . (3.9) 
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Obviously, they form the algebra 

= iF^, , [d^, u^] = 6; . (3.10) 

Next, the parallel displacement operator !P(x, x') of sections of the vector bun- 
dle S along the geodesic from the point x' to the point x is defined as the solution 
of the equation 

u"V^P{x,x') = Q , (3.11) 

with the initial condition 

r{x,x) = l, (3.12) 

where I is the identity endomorphism of the bundle S. Finally, we define the 
two-point quantity 

s^, = r~'v,r. (3.13) 

We remind, here, that we consider the case of a covariantly constant electro- 
magnetic field, i.e. 

V^F„^ = 0. (3.14) 
In this case we find it useful to decompose the quantity £^ as 

= -^^i^^ai^" + ^ • (3.15) 

By using this machinery we can rewrite the heat kernel as follows. First of all, 
the heat kernel can be presented in the form 

U{t\x, x') = exp i-t^) Pix, x')6{x, x') , (3.16) 

which can also be written as 

Uit\x, x') = rix, x')AHx, x') exp{-t^)6(u) , (3.17) 

where 6(u) is the usual delta-function in the normal coordinates (recall that u'^ 
depends on x and x' and u = when x = x') and ^ is an operator defined by 

£' = P~\x,x')A-kx,x')^Akx,x')r{x,x') . (3.18) 
As is shown in |l2l [TT|| the operator ^ can be written in the form 

= + (Vv + + ^v) 

= - (D^ + 4 - [Dy + Jy + , (3.19) 
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where = V^^ . 

Now, by using these equations and by recalling the formula in (13.151 ). one can 
rewrite the operator in (13.191) in another way as follows 

^ = - [xf^'D^Dy + Y^Df, + Z) , (3.20) 

where 

= rf,r]''\ (3.21) 

y = (V^Xf") + IX^^si^, , (3.22) 

z = ^;x^W,-^,x^^^, + (V^xn< + (V'')^v 

+ X^"V^< + X^^V/v • (3.23) 



4 Perturbation Theory 

Our goal is now to develop the perturbation theory for the heat kernel. We need 
to identify a small expansion parameter s in which the perturbation theory will be 
organized as e ^ 0. First of all, we assume that t is small, more precisely, we 
require t ~ s^. Also, since we will work close to the diagonal, that is, x is close to 
x', we require that u'^ ~ s. This will also mean that V ~ e ' and d, ~ s'^. Finally, 
we assume that F is large, that is, of order F ~ s'^. To summarize, 

u^'-s, F s'^ . (4.1) 
4.1 Covariant Taylor Expansion 

The Taylor expansions of the quantities introduced above have the form (up to the 
fifth order) [i2"l[TT1l 



1 

15' 



+ -^R^aApR^^u^u^u^u^ + 0{u^) , (4.2) 



+ ^R.aypR'ysWl^Uyu' + Oiu'), (4.3) 
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+ ^RapRysu^uPu-'u' + :^R^,aypR%\u"iJ'u-'u' + 0{u^) . (4.4) 
zoo JoU 

+ ^R^,ayp'R\u"j/u'' - \la'^p%^yU"l£u' 
Z4 o 

- ^R^aypR\^^F\u"^^u^u'u' + 0(u') . (4.5) 

We would like to stress that all coefficients of such expansions are evaluated at the 
point x'. Also note that, the expansion for is valid in the case of a covariantly 
constant electromagnetic field. 

4.2 Perturbation Theory for the Operator =Sf 

Now, we expand the operator ^ in a formal power series in s (recall that D ~ 
and M ~ e) to obtain 

oo 

(4.6) 

k=0 

where are operators of order £^~^. In particular, 

=2o = (4.7) 

^1=0, (4.8) 

= XfD^£), + y^D^+Z,, k>2. (4.9) 

where 

D' = g^'D^Dy, (4.10) 

and X^", and are some tensor-valued polynomials in normal coordinates u'^. 

Note that X'^^ are homogeneous polynomials in normal coordinates u'' and F 
of order s''. Similarly, ~ and ~ s''~^. Of course, here the terms Fuu are 
counted as of order zero. That is, they have the form 



C = (4-11) 
Zk = •P(4), k-2 + FapP^^)^ ^{u) + F,pFp^Pl^^l^^{u), (4.13) 
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where P^j)^ j,(u) are homogeneous tensor valued polynomials of degree k. 

By using the covariant Taylor expansions in (I4.2D . (14.51 ) and (14.31 ) we find the 
explicit expression of the coefficients 





— ' 


(4.14) 






(4.15) 




= H2 ajjU^U^ + L2 , 


(4.16) 






(4.17) 


^3 


— r/'^T//3 


(4.18) 


Z3 




(4.19) 




— C ^ Q^y^ U 14, , 


(4.20) 




= £4 al3yU"u^U^ + ajiySeU" u'' , 


(4.21) 


Z4 


— f^jji/'uP + Qj^y/^u" u7 uf' + O4 Qjjyii^Ku'^ uP u'^ , 


(4.22) 



where 







/ _ 




r^^ - 


—^R''{aliiFy)v 




1 

-^^RfiiJF p) , 


L2 = 


1 



1 1 2 

//3. = ^V^^R^ - iv^i? , (4.24) 
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G^4 Q-/3r(5e = '^R'^iaWlpR^r'^siF \A\e) > 

+ ^R,AyaR'''p + ^ARa^ + ^V.V^i? , 
L4a^yS = -■gQRfii»''/iR^7^F\y\s) - —R;i{a\A\/iR'^y''''iF\v\S) " 2^f^f^aR'^ / yiF\y\s) , 
O^aBydeK = T^TpRiiia' bR'^ y''' 6iF\v\eiF\A\K) ■ (4.25) 

J lb 

Here and everywhere below the parenthesis denote the complete symmetrization 
over all indices enclosed; the vertical lines indicate the indices excluded from the 
symmetrization. 



4.3 Perturbation Theory for the Heat Semigroup 

Now, by using the perturbative expansion (14.61) of the operator ^ and recalling 
that & ~ and t ~ e^, we see that the operator tD^ is of zero order and the 
operator t^k, k>2,h of (higher) order e^. Therefore, we can consider the terms 
t^k with > 2 as a perturbation. 

In order to evaluate the heat semigroup we utilize the Volterra series for the 
exponent of two non-commuting operators. Let X be an operator and 7 be a 
perturbation (say of order one in a small parameter). Then 

exp(X + y) = TexpX, (4.26) 

where 

^ 1 Tt Tl 

T = I + Yj J J dTk-i ■ Jdr, y(Ti)F(T2) ■ • ■ Y(Tk) (4.27) 

''"^ 

and 

y(r) = e'^Ye-^^ . (4.28) 
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By using the above series for the operator in (14.61 ) we obtain 

expC-?^) = T{t) expitD''-) , (4.29) 
where T(t) is an operator defined by a formal perturbative expansion 

J]r,(o, (4.30) 

k=0 

with Tk(t) being of order Explicitly, up to terms of fifth order we obtain 

(4.31) 
(4.32) 

V2(tTi), (4.33) 
V3(tT,) (4.34) 



Toit) 


= /, 




= 0, 


T2{t) 







-'! 






-'i 





1 T2 

.2 



J dT2 j dr, V2{tT,)V2{tT2) , (4.35) 




and 



Vkis) = e'^ ^te-'^^ . (4.36) 



4.4 Perturbation Theory for the Heat Kernel 

As we already mentioned above the heat kernel can be computed from the heat 
semigroup by using the equation (13.171) . By using the heat semigroup expansion 
from the previous section we now obtain the heat kernel in the form 

oo 

U(t\x,x') ~ r{x,x')A^'^{x,x')UQ(t\x,x')J^t''^^h(t\x,x') (4.37) 

where 

Uoitlx, x') = exp(tD^)6(u) , (4.38) 



/. G. Avramidi and G. Fucci : Non-perturbative Heat Kernel Asymptotics 16 



and 

bMx, x') = t-^I^U-Q\t\x, x')m)Uo{t\x, x') . (4.39) 

Thus, the calculation of the heat kernel coefficients reduces to the evaluation of 
the zero-order heat kernel Uo{t\x, x') and to the action of the differential operators 
Tk{f) on it. 

The zero order heat kernel UQ{t\x, x') can be evaluated by using the algebraic 
method developed in BUHl. First, the heat semigroup exp(?£)^) can be represented 
as an average over the (nilpotent) Lie group (13.101) with a Gaussian measure 

exp(?£)2) = (4;rO""^V(0 ^ ^^it exp |-iFM^,(OF + F£)^| . (4.40) 

R" 

where 

/ tiF 

■7(0 = det . (4.41) 

\ smh(nF) / 

and M{t) is a symmetric matrix defined by 

M{t) = iF coi\i{tiF) . (4.42) 

We would like to stress, at this point, that here and everywhere below all the 
functions of the 2-form F are analytic and should be understood in terms of a 
power series in F . 

Then by using the relation 

QX^{y'D^)6{u) = 6(u + k) , (4.43) 

one obtains 

UQ{t\x, x') = {Antyl^Jit) exp |-^m^M^^(Om''| , (4.44) 

which is nothing but the Schwinger kernel for an electromagnetic field on R." [2T1|. 

To obtain the asymptotic expansion of the heat kernel diagonal we just need 
to set X = x' {ox u = 0). At this point, we notice the following interesting fact. 
The operators fX/t, tVk{tT) and Tp,{t) are differential operators with homogeneous 
polynomial coefficients (in u^) of order . Recall that u ^ s, t ^ and F ~ s~^, 
so that tF and Fuu are counted as of order zero. Since the zero order heat kernel 
Uq is Gaussian, then the off-diagonal coefficients bk{t\x, x') are polynomials in u. 
The point we want to make now is the following. 
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Lemma 1. The off-diagonal odd-order coefficients b2k+i are odd order polynomi- 
als in u^, that is, they satisfy 

b2k+\{t\x,x) = -b2M{t\x,x), (4.45) 

andy thereforCy vanish on the diagonal 

4+iW = 0- (4-46) 

Proof. We discuss the transformation properties of various quantities under the 
reflection of the coordinates, u -u. First, we note that the operator D changes 
sign, and, therefore, the operator Xo = is invariant. Next, from the general 
form of the operator X/t discussed above we see that X/t ^ {-if-Lk- Therefore, 
the same is true for the operator V^itT), that is, Vk i-> {-l)'^Vk. 
Now, the operator T^it) has the following general form 

n = ■ ■ ■ Z ^'n,iVnitTi) ■ ■ ■ Vj„{tTj , (4.47) 

m=l ^ Lil=*: 

where the summation goes over multiindex j = (ji, . . . , j^) of integers ji, . . . , jm > 
2 such that |j| = j'l + ■ ■ ■ + j,,, = k, and C^j are some numerical coefficients. 
Therefore, the operator Tt transforms as Tt ^ {-VfTk. 

Since the zero-order heat kernel t/o is invariant under the reflection of coor- 
dinates u i-> -u, we finally find that the coeflicients bk transform according to 
bk ^ {-Vfbk- Thus, b2k are even polynomials and b2k+\ are odd-order polynomi- 
als. 

□ 

By using this lemma and by setting ;c = we obtain the asymptotic expansion 
of the heat kernel diagonal 

oo 

U^^\t) ~ {Antr"'^J{t) Yj t'^bl'^^it) , (4.48) 

k=0 

where the function J(t) is defined in (|4.41l) . Thus, we obtain 



(4.49) 
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4.5 Algebraic Framework 

As we have shown above the evaluation of the heat semigroup is reduced to the 
calculation of the operators Vk{s) defined by (14.361) . which reduces, in turn, to the 
computation of general expressions 

e^^^'M^' • • -M^-D^, ■ ■■'D,„e-'''" = Z^'(5) ■ • ■Z^"(5)A^,(5) • • -A^S^) , (4.50) 

where 

Z\s) = e'^'u'e-'^' . (4.51) 
A,{s) = e^'^'n.e-''^ . (4.52) 
Obviously, the operators and Zy form the algebra 

[A,is),Z\s)] = 6;, [A,{s),A,{s)]=iF^y, [Z^(^), Z^(^)] = . (4.53) 

The operators A^(s) and Z^is) can be computed as follows. First, we notice 
that Af,(s) satisfies the differential equation 

dA^is) = M^^A^is) , (4.54) 

with the initial condition 

Hereafter Ad^i is an operator acting as a commutator, that is, 

Ad^A.is) = [D\A^is)] . (4.55) 
The solution of eq. (14.541) is 

A^(5) = expis Ad^2)D^ , (4.56) 
which can be written in terms of series as 

Now, by using the algebra (13.101) we first obtain the commutator 

[D\DJ = -2iF.M\ (4.58) 
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and then, by induction, 

{M^.f = {-2ifF,,,F"',, ■ ■ ■ = [i-2iF)%,n" . (4.59) 
By substituting this result in the series (14.571) we finally find that 

A,is) = , (4.60) 

where 

*P(^) = exp{-2siF) . (4.61) 
Similarly, for the operators Z^(s) we find 

Z'^is) = exp(5Ad2)2)M^ = y 77 (Ad2)2f u" . (4.62) 

k=Q 

Now, by using the commutators in (13.101) . we find 

Ado2U^' = [d\ m^] = 2ir , (4.63) 
and then, by induction, we obtain, for k >2, 

iAd2j2f u^" = 2[i-2iFf-^r''Da . (4.64) 
Thus the operator Z'^is) in (|4.52l) takes the form 

t 
kl 

This series can be easily summed up to give 

Z'^(5) = u'' + a>'"(s)Da , (4.66) 

where 

l-expi-2siF) sinh(^zF) 

ais) = — = 2expi-siF) — . (4.67) 

iF iF 

Now, by using (|4.61l) and (14.671) we obtain 

Q.-\s) = ^iF [coth(5/F) + 1] = ^ [M(s) + iF] . (4.68) 



Z^(s) = - 2sir + 2 y ^[(-2iFf-^r"Da . (4.65) 

k=2 
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We will need the symmetric and the antisymmetric parts of Q. ^{s). By recalling 
that the matrix F is anti-symmetric it is easy to show 

= \m,M . (4.69) 

"i(^) = ^'^/^v , (4.70) 

Here and everywhere below the square brackets denote the complete antisym- 
metrization over all indices included. 

For the future reference we also notice that 

Qr\s)QF{s) = ^-\s) = exp(2siF) , (4.71) 

Finally, we define another function 

cD(^) = '¥(s)Q-\s) = (Q~\s)f = ^ [M(s) - iF] . (4.72) 
It is useful to remember that the functions FQ. and Of! are dimensionless. 



4.6 Flat Connection 

Next, we transform the operators to define new (time-dependent) derivative 
operators by 

D,{s) = n;l{s)Z\s). (4.73) 
By using the explicit form of the operators and we have 

= V^ + ^M^p(s)u'' . (4.74) 

Since the operators commute, the operators D^,(s) obviously commute as 
well. In other words the connection D^, is flat. Therefore, it can also be written as 

D^is) = e-®^'^Vf,e®^'\ (4.75) 

where. 



0(^) = -u^M^yis)u'' (4.76) 



1 
4 

Now, we can rewrite the operators A^(5) and Z^(s) in (14.601) and (14.661) in terms 
of the operators Dfj{s) 

Z^{s) = Q!"'{s)Da{s) . (4.77) 
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5 Evaluation of the Operator T 

The perturbative expansion of the operator T is given by the eq. (14.301) . with 
the operators being integrals of the operators Vk{s) and their product. Thus, 
according to (|4.33D - (|4.35I) . to compute the operator T up to the fourth order we 
need to compute the operators V2{s), V-i{s), V^^is) and V2{si)V2{s2). 



5.1 Second Order 



Now, by using the explicit expression for ^2 given by eqs. (14.91 ). (14.161) and (14.231) . 
utilizing the results of the Section 3, exploiting eqs. (14.771 ). (IB. 21 ) and (IB. 31 ). using 
eqs. (14.611 ). (14.671) . (14.711 ) and (14.721 ) after some straightforward but cumbersome 
calculations we obtain 



where 



1 

6^ 



A^(2) = 



jY6 

(2) 



cr6y 
(2) 



1 



fIV 



fyPcrSy 



'(2) 



(5.1) 



(5.2) 



(5.3) 
(5.4) 
(5.5) 



Note that all these coefficients as well as the operators depend on the time 
variable s. We will indicate explicitly the dependence of various quantities on the 
time parameter only in the cases when it causes confusion, in particular, when 
there are two time parameters. 



5.2 Third Order 



Similarly, by using the explicit expression for ^1, given by (14.91 ). (14.191) and (14.241) . 
utilizing the results of the Section 3, exploiting eqs. (14.771 ). (IB. 21 ) and (IB. 31 ). using 
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eqs. (14.611 ). (14.671) . (14.711 ) and (14.721 ) after some straightforward but cumbersome 
calculations we obtain 



where 



^(3) 



(5.6) 

(5.7) 
(5.8) 



(5.9) 
(5.10) 
(5.11) 



Here again, for simplicity, we omitted the dependence of the coefficient functions 
and the derivatives on the time variable s. 



5.3 Fourth Order 
5.3.1 Operator V^is) 

By taking into account the definition of ^4 in (14.91) by using eqs. (I4.20I) - (I4.22I) . 
(14.771) . (IB. 21) and (IB. 31) . and the explicit form of the functions *F and Q, we obtain 

+ Y'^P;''D^DpD,D,D, + S 1^;"''D^DpDfi,D,D;, , (5.12) 
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where 



(5.14) 



300 

X 



, 1 



^ 240 



^ 240 



(5.15) 



- y^i^/^^^'^^^-'Qa^-fi/fl/i^/ [3(5/) + 13^/)] ^^u^ 



(5.16) 



20 

125^^'^^/) + 55^'^^/)] . (5. 17) 
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5.3.2 Operator V2{si)V2{s2) 

Next, we need to compute the product of two operators V2{s) depending on dif- 
ferent times s\ and S2 by using the eq. (I5.1I ). To simplify the notation we denote 
the derivatives D^{sk) depending on different times Sk simply by D®. To present 
the product V2{s\)V2{s2) in the "normal" form we need to move all derivative op- 
erators D^'^ to the right and all coordinates to the left. In order to perform this 
task we need the commutator of the derivative operator D^^^ with the coefficients 
of the operator V2{s2)- First, by using the commutators listed in Appendix B we 
obtain the relevant commutators 

nf\,Xs2)D'il ■ ■ ■ D'^l , (5.18) 

«(«-l)^'W(^2)Z)i^••<^) 
+^/^'V.(^2)Z)i^■■Z);^^_^), (5.19) 
np^^\,,{s2)D^^---D%, (5.20) 

where 

A = {n''is-\R''}^d''%A, (5.21) 

h'^^A = ^R^(a''fi)^"'^'%.^rAU\ (5.22) 

p"''a = -^i?^^/^f2'^('Q^I''[C^ + 7¥/)]o,, . (5.23) 

Next, by using the expression for the operator V2{s) in (15.11) and the non- 
vanishing commutators in (|5.18I )- (I5.19I ) we obtain 

^2(^1)^2(^2) = ^R^ + ^R[V2iSi) + V2iS2)]+LiSuS2), (5.24) 

where 

L(.i,.2) = XZc™--(.,,.2)/)i?---^^^^^^^^ , (5.25) 

k=\ n=Q 



Pi?---^iiUb(^2)] = 

H!---DZP?2M = 
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and 



Cl',^ = 2A^('^2)('Vl)<2)('^2) + 2P'(^)(.i)r,(.^2), 

Cjl) = 2Pg(.iWf2)(^2) + 3H^pf(5i)r.(^2), 

CIS) = 2W(2)''(^i)A^f2)(^2) + 4!2jr(^i)rt^2), 

C(t[f = 2Q^^;\sON^2)(^2) , (5.26) 



CTi) = 2A^;^)(5i)Pg;(52) + 2Pg(5i)/i^V^2) + 6Wg^(^i)^^V(^2), 

CZ) = 2Pj^(.0/t2';(^2) + 3Wg^(50/^^%(52) + 12!2j/'(^i)g^V^2), 

= 2W;g^(5i)PS(^2) + 4!2j/'(5i)/i^'^.(.2), 

^jrjp. ^ 2eg^^(.i)i^2^(^2) , (5.27) 



C;^';37 = 2Ar;^2)(5l)Wj7(52) + 2/>g;(5i)K%(52), 

^jpc. ^ 2P;5(5i)Wj7(52) + 3Wjf (5i)p^%(52) , 

cj^p = 2w;f(5i)wg7(52)+4(2g;^^(5i)p^%(52), 

^«<.. = 2!2g^^(.0<27(^2) , (5.28) 



c™ = A^a)(^i)e?27''(^2) , 

cSr = <)(^i)e^27^(^2) , 

C|JrP-^ = H^(t(50e^27^(.2) , 

^a,ySpa., ^ Q"^y\s,)QP^;'{S2) . (5.29) 



6 Generalized Hermite Polynomials 

Thus, we reduced the calculation of the asymptotic expansion of the heat kernel 
to the calculation of the derivatives Di^{s) of the zero order heat kernel UQ{t\x, x') 
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given by (14.441) . The needed derivatives of the zero order heat kernel can be ex- 
pressed in terms of the following symmetric tensors 

9{,,...,Ss) = U^\t\x,x')D^,(s) ■ ■■D^Ss)Uo(t\x,x') (6.1) 

and 

S,,..,,„,,..^„(5i,^2) = U^\t\x,x')Di';> ■■■D^;^D^^^ ■■■Df%(t\x,x') , (6.2) 

where we denoted as before D|f^ = D^(sk). 

We recall that the derivatives O'j}-' and D^y^ do not commute! Also, Uq is a 
scalar function that depends on x and x' only through the normal coordinates m^. 
The derivative operator Z)^(5) is defined by (14.741) . and, when acting on a scalar 
function is equal to 

d 1 

= (6.3) 

where the tensor M^y{s) is defined by (14.421) and the function Q(s) is a quadratic 
form defined by (14.761) . 

Therefore, by using the explicit form of the zero order heat kernel (14.441) we 
see that the tensors 'Kpi.../j„(^) can be written in the form 

'H^,..^„(5) = exp{0(O - 0(5)}—- • • • — - exp{0(5) - 0(0} , (6.4) 

The tensors 'H^^.-.^^X^) are polynomials in u''. They difl'er from the usual Her- 
mite polynomials of several variables (see, for example, [16]) by some normaliza- 
tion. That is why, we call them just Hermite polynomials. The generating function 
for Hermite polynomials 

oo 

s) = y -e' ■ ■ ■ e-"H,,...„xs) (6.5) 

can be computed as follows 

•K(^,5) = exp {0(0- 0(5)1 exp|^^^jexp {0(5) -0(0}, 

= exp|irAaM5)[f^ + 2i<'^]| , (6.6) 
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where 



A(^) = \[M{s)-M{t)\ 

(t - xMF^ 

(6.7) 



2l 

1 iF sinh[(? - s)iF^ 

2 sinh(?/F) sinh(sjF) 



By expanding the exponent in ^ we obtain the Hermite polynomials explicitly. 
They can be read off from the expression 

[f] 

• • • = g ^(;^) {^^K,{s)e)' [e^,.is)u'^r' . (6.8) 

For convenience some low-order Hermite polynomials are given explicitly in ten- 

sorial form in Appendix A. 

Similarly, the tensors 'E.y^...v^m...^J^s\, S2) can be written in the form 

^Vv-V,nUV-Hni^l^ ^2) = BXp [0(0 - ©(^j)] (6.9) 

• • TIT exp [©(^i) - 0(^2)] [®(*2) - 0(0] 

They are obviously polynomial in as well. We call them Hermite polynomials 
of second kind. The generating function for these polynomials is defined by 

oo - 

S(^, 77, Su S2) = y — -r • • • • • • ?/'"Sv,..v„^,..^„(5l, 52) , (6.10) 

mini 

m,n=0 

and can be computed as follows 

S(^, 77, 51, 52) = exp {0(0 - 0(5i)} exp l^e-^^ exp {©(sO - 0(52)} 

X exp l?;^^ j exp {©(52) - 0(0) , (6. 1 1) 

= exp |i^'^A„^(5i)(^ + 2u^) + ^rfA,y(s2)(ri' + 2u^) + fA,M2)r]" 
Notice that 

5(^,77,51,52) = 'H(^,5i)<K(77, 52)exp{^''Ap<,(52)/7") (6.12) 
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This enables one to express all Hermite polynomials of second kind S(„)(5i, S2) 
in terms of the Hermite polynomials "K(„,)(5i), ^(i){s2), and the matrix A(52). 
Namely, they can be read off from the expression 



1 Off-diagonal Coefficients bk 

By using the machinery developed above, we can now write the coefficients of the 
asymptotic expansion of the heat kernel in terms of generalized Hermite polyno- 
mials. We define the following quantity 

1 

b2Xl){t\x, X') = J dT[N[^^itT)9{AtT) + PllpT)nHys{tT) + W'g^\tT)^^ys{tT) 


+ &^2^\tmp^ys{tT)\ . (7.1) 

Then, by referring to the formulas (15.11 ), (15.61 ). (15.121 ) and (15.241 ) and by using the 
following formula for multiple integrals 

h T„ T2 b 

J dTn J dTn-l ■ J dTj{T,) = J dz (b - t)"-' /(t) , (7.2) 



a a 



we obtain 



b2{t\x,x') = lR + b2xi)it\x,x'), (7.3) 



1 



b3it\x, x') = r'l^ J dr[N[y^itT)'HAtT) + PllpT)nHys{tT) + W'^^\tT)^^ys{tr) 







+ Q';)'\tT)nHp^ysitT) + Yl^;''\tT)Up<ry6(tT)] . (7.4) 



/. G. Avramidi and G. Fucci : Non-perturbative Heat Kernel Asymptotics 29 



bA{t\x,x') = + lRb2xi)(t\x,x') 



1 



4 4 1 -^2 

YjYj fdr2 f jTic;;;,f'^'--^*ari,?T2)s^,..p„,,...,aTi,?T2). 



+ 



+ 



(7.5) 



8 Diagonal Coefficients bk 

In order to obtain the diagonal values bf'^^it) of the coefficients bkit\x, x') we just 
need to set m = in eqs. (17.31) . (17.41) and (17.51) . For the rest of this section we will 
employ the usual convention of denoting the coincidence limit by square brackets, 
that is, 

[/(«)]'"'= /(O). (8.1) 

By inspection of the equation defining the generalized Hermite polynomials in 
Appendix A one can easily notice that, in the coincidence limit, all the ones with 
an odd number of indices vanish identically, namely 

n.-P..,] =0. (8.2) 

By using the last remark we have the following expression for the coincidence 
limit of (1731) . i.e. 

bf''it) = ^R + bi%(t), (8.3) 

where 

1 

<a)(0 = J rfT[/'g')(?T)'K,,(rT) + e^^'''(rT)'Hp^,,(fT)]'"' . (8.4) 



By using the explicit form of the coefficients ^(2), 2(2) and the generalized Hermite 
polynomials in Appendix A, we obtain 

^Zfp) = Jm^^vm^ap + Jiirm^y + Ji^rm,, , (8.5) 
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where 





+^ (5/ + 3¥/) a«^^/-¥/A^,A,,)| , (8.6) 
1 





1 



Ji^fit) = J dTQ.^y'i"'^^A''s ■ (8.8) 



Here all functions in the integrals depend on tr. 
Next, we introduce the following matrices 

car \ ^ 1 exp[(r - 2s)iF] - exp(-tiF) 

Jl(s) = a(s)A(s) = ■ . , . (8.9) 

2 sinh(rjF) 

^ coth(?/F) cosh[(? - 2s)iF] 

iF iF siian{tiF) 

T(s) = Q-\s)-^^(s)A(s)-^A(s) 

= - I B/F coth(rjF) + — — — cosh[(? - 2s)iF]\ . (8.1 1) 
8 \ sma(tiF) ) 

Then, by using the relation 

Q.{s)A{syi'{sf = a^(s)A(s) = J(^(s) (8.12) 
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we obtain 



1 



+\ {^,^"{tTW\{tT) + 3^(/{tT):R,f{tT)) I , (8.13) 
1 

W« = \^dT:R^>'\tT)=^-^\ (8.14) 



WW = - j6?T:?l^%T) = -^|^-coth(f/F)J^ . (8.15) 



Unfortunately the integral J{i)"^^v can not be computed explicitly, in general. 

As we already mentioned above all odd order coefficients b2k+\ have zero di- 
agonal values. We see this directly for the coefficient hi,, which is given by (17. 4D . 
That is, by recalling the formulas in (15.71) through (15.1 II) and the remark (18.21) we 
have 

bf\t) = Q. (8.16) 

Finally, we evaluate the diagonal values of fourth order coefficient hi, given by 
(17.51 ). It can be written as follows 

bT\t) = ^R' + ^RKZit) + K2,it) + <J)(0 . (8.17) 
By noticing that for odd n + k, the diagonal values of the coefficients C(^n,k) vanish, 

= 0, (8.18) 

and by using the explicit form of Hermite polynomials and the generating function 
(16.121 ) we obtain 

1 

= J jT{i"(^)(?T)A,,(?T) + 3[e;^f(fT)]'"'A(,,(?T)A,,)(?T) 



+155;f' ^(rT)A(,,A,,(?r)A,^)(rT)} , (8.19) 
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(2) 



+ 
+ 
+ 





+ 12e^^^^(rO^<--V.(r,)) (Ai^A?) + 2A2)Ag;) 
+ 6[Pj«(rO]%<-'^^(T,)ASA^ 

2<^(r,)[/-;(r.)]^'^^ (3A;>;>2) + 12ASa(?aS) 
4<2?r(T0p<-'^^(r,) (9ASaSa^S + 6AgAgA^)) 

err (roe?r(-^)(9ASASA?x? + 

+ 24Ag)A>^^)Aj)} , (8.20) 

where the superscript on the matrix A denotes its dependence on either tTi or tT2. 
We see that the scalar curvature appears only in the term b^^^-p). Now, the 

term ^4 (2)(0 only contains derivatives of the curvature and quantities which are 
quadratic in the curvature with some of their indices contracted. It has the follow- 
ing form 

+ A5„^(0V«V^i? + A^^^COV^V^^R^^ + ^5,^(0V-V^7?^^ , (8.21) 
Here the tensors a^'\s) are functions that only depend on F (but not on the Rie- 
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mann curvature) defined by 

3 13 

Aayp{s) = —^(ay^fiU + -^^{ay^\A\/3) , (8.22) 

a^'^aV<5(^) = ^S„osS^^(31(m)'^^> + 65A'^^)-lM^''S,o8S^,) 
187 3 1 

+ ^%,^^^5)^'^)„ + ^(S«o8^<^^,/^ + %r^5)^'^/^) ' (8-23) 
3 1 

cP\tjfir(s) = -^^(ap^)n + gQ^CajS^^ir) > (8.24) 

d'^'^'' aixPy{s) = --S{aji^y)n ~ g^(ffyS^|A<|r) , (8.25) 

+ ^^ipy^^^s^^'K, (8.26) 
9 3 

(i^^^aPnv{s) = 2Q^{aP'^\M\v) ~ ~^^n(a^M ' (8.27) 

3 

a^''\pnv{s) = --S^aji^W) • (8.28) 

All functions here are evaluated at the time s (unless specified otherwise). 

The term ^4 1'3)(0 only contains quantities which are quadratic in the curvature 
with none of their indices contracted. It has the form 

+ ^%pp^t)f^''R"''"' + ^ZfiitW^R"" + ^ZfiitW'R"'' , (8.29) 

where D^'j...^„(0 are some tensor-valued functions that depend on tF. They have 
the form 

D2...,„(0 = JdT2 j dT,dl:l..^StTutT2). (8.30) 


To describe our results for the tensors d^'^^ we define new tensors 

= V + /^V' (8-31) 

- ^Q/''n,^£(i)/>0,,A^A;,e + ^ilp^'ila''i'a'S^s)/^A,,A,,A,^ (8.32) 
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- ^(A,,Ar,^ytTi){!BsaMypytT2) 

- K,(tTi)ijls(^Jli^ir,)MypytT2) 

- \ \^Kr,{tT^)(K^A<^r){tT2) + ^ A^,(rTi)(A;,^A^^)(?T2)| 

X (^l/^Q,-fi(l)/>0^,)(^T2) 

+ ^|(At/.Aw)(^'^i)(A6TA^^)(?T2) + 8A,^(?Ti)(A,;,A^^A^a)(?T2) 

+ ^(A,,A,,A,,A;,,)(?T2)|(i^/^r2,^^/£(3)p-'))(?T2) . (8.33) 
Then the tensors d^'^^ have the form 

+ ^(S^vM„J(?Ti)(S,,Q^^y(?T2) 

+ i(Q/'Q/)M,J(rri)(j?l(,|,|J?l<,),M;;)(?T2) 

+ ^(f2(^':?lv),M„^)(fTi)(Q,(''D/£(i)/0^,A.;,)(?T2) 

- ^{(Sa^Myp)(?Ti)(A,^A^^)(?T2) 

+ 4(Q/-t2/)M^p)(?Ti)(A,Aa.A,^)(rT2))(il/'f2/^„''£(3)/)(rT2) 

+ ^(f2/'Qv'^^a''fi(3)/^)(?Tl)%,p<.,^;,(?Ti,?T2), (8.34) 
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- ^(S^crM„p)(?Ti)j?l(^,)(?T2) 

- \^iMy)^tTi)[Sp^Map){tT2) 

+ ^(^(/^v)A(A;^)(^Tl) + ^A,^(rTl)(:?lv(;,A;^,)e£(7)p')(/'T2) 

+ 4Qv^'£(7)/^aTl)<S«^p,,,(?T2) , (8.35) 

Jo 

d'tU(tTutT2) = -3^(n/'£(7)/))(rTi)(0<,,:?l^,)(rT2)+ ^:?l(;,v)(?Tl):?l(<,^)(rT2) 

+ ^^lv^'£(7)/^(?Ti)(j?l^,A,,£(7)„'^)(?T2) . (8.36) 
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+ i(Q/^Y/>)(fTi)(Q/'f2/£(i)p''>0,,A,,A,,)(?T2) 

+ i(D/Q/¥,''£(3)p^))(fTi){2A,,(rTi)(j?l,,$^,)(rT2) 

- (j?l^vA..A,e)(?Tl) - 4A,,(?Ti)(%^^v,A,^)(?T2)) 

- ^{(j?l;.vA,,A,,)(rT2) + 4(n,('^^/))(rTi)(A,A^A,.)(?T2)) 

X (Q/'Q/^,''6(3)p^^)aT2) . (8.37) 

+ ^(a,(''¥;^))(?Ti)(o,,j?[^,)(?r2) 

2 2 

- ^(^l/'fi(7)a'^^)(^Tl)(%^:?lvA.)(^T2) 

- ^(a^'V)aTi)(Q(''/^A,,A,,)aT2) . (8.38) 

<l/^^l'^^2) = -2(Q,(^^/))(?Ti)(0<^J?l^,)(?T2) + 2j?l^v(?Ti)j?l,^(?T2) 

+ 4(a/^¥/))(rTi)(^„'-:?l;j,A,,)(rT2) . (8.39) 

9 Conclusions 



In this paper we studied the heat kernel expansion for a Laplace operator acting 
on sections of a complex vector bundle over a smooth compact Riemannian man- 
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if old without boundary. We assumed that the curvature F of the U{\) part of the 
total connection (the electromagnetic field) is covariantly constant and large, so 
that ?F ~ 1, that is, F is of order ? In this situation the standard asymptotic 
expansion of the heat kernel as ? — > does not apply since the electromagnetic 
field can not be treated as a perturbation. 

In order to calculate the heat kernel asymptotic expansion we use an algebraic 
approach in which the nilpotent algebra of the operators plays a major role. 
In this approach the calculation of the asymptotic expansion of the heat kernel is 
reduced to the calculation of the asymptotic expansion of the heat semigroup and, 
then, to the action of differential operators on the zero-order heat kernel. Since 
the zero-order heat kernel has the Gaussian form the heat kernel asymptotics are 
expressed in terms of generalized Hermite polynomials. 

The main result of this work is establishing the existence of a new non-pertur- 
bative asymptotic expansion of the heat kernel and the explicit calculation of the 
first three coefficients of this expansion (both off -diagonal and the diagonal ones). 
As far as we know, such an asymptotic expansion and the explicit form of these 
modified heat kernel coefficients are new. 

We presented our result as explicitly as possible. Unfortunately, some of the 
integrals of the tensor- valued functions cannot be evaluated explicitly in full gen- 
erality. They can be evaluated, in principle, by using the spectral decomposition 
of the two-form F, 

[n/2] [n/2] 

F = J]B,Ek, F^ = -J]bIu,, (9.1) 

k=l k=l 

where are the eigenvalues, are the (2-dimensional) eigen-two-forms, and 
Ilk = -El are the corresponding eigen-projections onto 2-dimensional eigenspaces. 
Then for any analytic function of tiF we have 

[n/2] . [n/2] 

fitiF) = fm)-(Uk + iEk) + 2 /(-r5,)-(n, - iEk) . (9.2) 

k=l k=l 

However, this seems impractical in general case in n dimensions. It would sim- 
plify substantially in the following cases: i) there is only one eigenvalue (one 
magnetic field) in a corresponding two-dimensional subspace, that is, F = BiE\ 
(which is essentially 2-dimensional), and ii) all eigenvalues are equal so that 
= -I (which is only possible in even dimensions). We plan to study this 
problem in a future work. 
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The work carried on in this paper can find useful appUcations in various fields 
of theoretical physics and mathematics. For instance, our results can be applied to 
the study of the heat kernel asymptotic expansion on Kahler manifolds. The com- 
plex structure on Kahler manifolds is a parallel antisymmetric two-tensor which 
plays the role of the covariantly constant electromagnetic field. This subject is 
also interesting, in particular, in connection with String Theory. 

Appendix A. Hermite Polynomials 

The Hermite polynomials are defined by 




(A.1) 



They can be computed explicitly as follows. First, let 



(A.2) 



and 




(A.3) 



Then 



'K(„)(^) = (A + 5)''-l. 



(A.4) 




(A.5) 



{B,A\ = C, [A,C] = [5,C] = 0. 



Lemma 2. There holds, 




(A.6) 
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Proof. Notice that is the generating functional for (A + B)". Now, by using 
the Baker-Hausdorff-Campbell formula 



expanding both sides in t and computing the Taylor coefficients of the right hand 
side we obtain the eq. (IA.6I ). □ 



By using this result we obtain an explicit expression for ( IA.4I ) 

n\ 

2^k\{n-2k)\ 



[f] 

k=Q 



By setting A = we immediately obtain the (diagonal) values of Hermite polyno- 
mials at a = 

[K.-,2,,..r' = 0, (A.8) 

r -idiag (2^)! 

|^'7l^,...jU2„J = ~^J^^^(Mifi2 ' ' ' ^M2n-if2n) ■ (A. 9) 

We list below a few low order Hermite polynomials needed for our calculation 

= u (A.10) 

-H^, = A^,,u\ (A. 11) 

'J^Mifi = A(^„^2) + Api^A^iX"^ , (A. 12) 

'J^fiiii2fi3 ~ 3A(p,yj2Ayj3)(j,M + A^jaAyjj^ApjyM i£u7 , (A. 13) 

'J^lllM2M3f4 - ^^(piM2^M3f4) + 3A(p,^jA^3|Q,|Ap^)^M l/ 

+ ApjQ,Ap2/jA^3yA^^5M u^uJu^ . (A. 14) 

'J^fiiii2tJ3M4M5 ~ 15A(p,p2A^3^4Apj)Q,M + 5A(pj^2Ap3|(i-|Ap4^g|Apj)ya u^u^ 

"I" A^jQ.Ap2/3Ap3yAp45Apj;yM "u^uJu^u^ , (A. 15) 

^flllJ2tJ3l^4MSM6 ~ ^5^(pifl2'^M3M4^I^Sl^6) ^^'^(MlM2^f^3fi4'^MsM^f6)/i^ 

+ 1 5 A(^ Ap3 \a\ ^M4 A/^s Irl ^1^6)6^"^^ 

+ A(p,|Q,|Ap2te|A/j3|y|Ap4|5|Ap5|;^|Apg),M"M^M''M''M''M' , (A. 16) 
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We list below some of the generalized Hermite polynomials of second kind. 
Now we have two sets of Hermite polynomials that depend on the quadratic forms 
A at two different times, ^'i and si- Let us define 

<K(„)(^i) = r ■■■f'"Hr-/.„(^i)> (A.17) 
•H(„)(52) = ^' • • • rf"H,,..,,Xs2) , (A.18) 

and 

Ais2)=eKM)^ ■ (A.19) 
Then from eq. (|6.13l) we obtain the quantities S(m,n) that we need in our calcula- 
tions 



S(0,l)(5l,52) = 




(A.20) 




A(*2) + 'H(i)(5i)'K(i)(52) , 


(A.21) 


S(2,l)(^l,^2) = 


2A(52Wi)(^i) + 'Ha)(s2)%2)(si) , 


(A.22) 


S(3,l)(5l,52) = 


3A(52)'K(2)(5l) + 'K(l)(52)'K(3)(5l) , 


(A.23) 


S(4,l)(*l,*2) = 


4A(*2)^(3)(^l) + 'H(i)(52)'H(4)(*l) . 


(A.24) 


S(0,2)(^l,^2) = 'H(2)(52) 




(A.25) 


S(l,2)(*l,*2) = 2A(52)'K(i)(52) + 'H(2)(52)'H(i)(5i) , 


(A.26) 



E(2,2)(*l, *2) = 2X\S2) + 4A(52)'K(i)(52)'H(i)(5i) + 'K(2)(^2)'K(2)(5l) (A.27) 
S(3,2)(*l,*2) = 6A2(52)'K(1)(5i) + 6A(52)'H(i)(52)'H(2)(5l) 

+ 'H(2)(52)'K(3)(5i) , (A.28) 

S(4,2)(*l, *2) = 12A2(52)'K(2)(5i) + 8A(^2)'?^(l)(^2)'?^(3)(5l) 

+ 'H(2)(^2)'H(4)(^i) . (A.29) 

S(0,3)(^l,^2) = 'K(3)(^2) (A.30) 
S(l,3)(*l,^2) = 3A(^2m(2)(52) + 'H(3)(52)'H(i)(5i) , (A.31) 
S(2,3)(*l,*2) = 6A2(52)'K(1)(52) + 6A(52)'K(2)(52)'H(i)(5i) 

+ 'H(3)(^2)'H(2)(*i) , (A.32) 

S(3,3)(*l,*2) = 6A3(^2)+18A2(^2)'H(1)(52)'H(1)(5i) 

+ 9A(*2)'K(2)(*2)'H(2)(5i) + 'K(3)(52)'H(3)(5l) , (A.33) 

S(4,3)(^l,*2) = 24A3(52)'H(1)(5i) + 36A2(^2)^(1)(^2)'H(2)(^i) 

+ 12A(52)'K(2)(^2)'H(3)(^l) + 'K(3)(^2m(4)(5l), (A.34) 
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S(0,4)(5l,52) = 'H(4)(52) (A.35) 

3(l,4)(5l,52) = 4A(52)'H(3)(52) + 'H(4)(52)'H(i)(5i) , (A.36) 

3(2,4)(5l,52) = 12A2(52)'H(2)(52) + 8A(52)'H(3)(52)'H(l)(5l) 

+ 'H(4)(52)'H(2)(5l) , (A.37) 

3(3,4)(5l,52) = 24A^(52)'H(l)(52) + 36A'(52)'H(2)(52)'H(l)(5l) 

+ 12A(52)'H(3)(52)'H(2)(^i) + 'K(4)(52)'H(3)(5i) , (A.38) 

S(4,4)(5l, S2) = 24A^(^2) + 96A3(52)'H(1)(^2)'H(1)(^i) + 72A2(52)'K(2)(52)'H(2)(5i) 
+ 16A(52)'K(3)(52)'K(3)(5l) + 'K(4)(52)'K(4)(5l), (A.39) 

The coincidence limit of the quantities „) ,with m + n odd, vanishes identi- 
cally 

[S(^,„)(5i , 52)]'"' = , if (m + n) is odd (A.40) 

By recalling the coincidence limits of the Hermite polynomials we obtain the 
following 



[S(ij)(5i,52)]'"^^ ■■ 


= A(S2), 




(A.41) 


[S(3,i)(5l,52)]'"^^ ■■ 


= 3A(5i)A(52), 




(A.42) 


[S(0,2)(5l,52)]'^^^ : 


= A(52), 




(A.43) 


[S(2,2)(*l,*2)]'^^^ ■■ 


= AisMs2) + 2A\s2), 




(A.44) 


[S(4,2)(5l,52)]'^^^ ■■ 


= 3A\si)A(s2) + l2A(s,)A\s2) , 




(A.45) 


[S(l,3)(5l,52)]'^^^ : 


= 3A\s2), 




(A.46) 


[S(3,3)(5l,52)]'*'^^ : 


= 9A(si)A\s2) + 6A\s2), 




(A.47) 


[S(2,4)(5l, S2)f^^ ■■ 


= 3A(sOA\s2) + 12A^(52) , 




(A.48) 


[S(4,4)(ji, 52)]*^^ ■■ 


= 9A\si)A\s2) + 72A(5i)A^(52) + 24A^(52) 


. (A.49) 



Appendix B. Commutators 

Lemma 3. Let and be operators satisfying the algebra 

[D^, m"] = 61 , [D^, Dy] = [M^ m"] = . (B.l) 

Then 

[D^r--D,„,u''] = n6P(,,D^,---D,^^ (B.2) 

+ 2nM'^5"VD^,---D^„). (B.3) 
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Proof. Let X{0 = ^''D^ and 

^{t) = up] = (e'^'^^^u^e-"^'^^ - u'') e"^^^^ , (B.4) 

Then 

00 I. 

e'^i^^u^e~''ii^ = J]UAd,,,,fu'^. (B.5) 

k=0 

By using the commutation relation in (IB.ll) we have 

[X(o,un = e (B.6) 

and, therefore, 

e'^^^^u^e-'^^^^ = uf + . (B.7) 

Thus 

i/(t) = t^e''^^^^ . (B.8) 
By expanding in Taylor series both sides of the last equation we obtain 

Z TTTTTT^" • ■ ■ K. • • • ^^-.)' "1 = Z TT^" ' ' ' ^''''^(m.D,._ ■ ■ ■ D,,^,, . 

k=0 ^ ^' k=o ^• 

(B.9) 

Now by equating the same powers of t in both series we obtain the claim (IB.2|) . 

The second relation can be proved in a similar manner. We introduce, in this 
case, the following generating function 

^-(f) = . (B.IO) 

By the same argument used in the proof of the first relation we obtain that 

^^^^(t) = [e'^*^^^ , uPu''] = It^Ve'^^^^ + t^e^" ■ (B . 1 1 ) 

Now, as before, by expanding the last equation in Taylor series and equating the 
same powers of t we obtain the claim (|B.3I) . □ 
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